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ABSTRACT 


Simultaneous  direct  measurements  taken  in  eleva- 
tion, azimuth  and  range  from  independent  instru- 
mentation sources,  when  aggregated  and  reduced 
can  often  yield  sensible  estimates.  However, 
strong  aberrations  and  inconsistencies  in  the 
estimates  derived  from  data  reduction  have  been 
known  to  occur.  Two  particular  factors  arc 
investigated  in  order  to  partially  explain  theo- 
retically these  elusive  deviations  and  to  estab- 
lish some  cause  and  effect  relationship.  The 
first  set  of  factors  relates  to  the  relative 
geometric  siting  of  the  instruments  as  well  as 
the  type  of  measurements  recorded  • Typical  mea- 
surements which  would  be  conducive  to  such 
aberrations  are  discussed.  The  second  set  of 
factors  relates  to  the  selection  of  a frame  of 
reference  for  azimuthal  measurement.  It  is  shown 
that  coordinate  selection  can  play  for  certain 
types  of  measurements  a significant  role  in 
multiplying  the  GDOP  by  severalfold  from  a ratio- 
nally established  minimum  basis. 


estimates  derived  from  data  reduction  have  been 
known  to  occur.  In  many  instances,  the  cause 
of  these  deviations  cannot  be  correlated  with 
such  factors  as  the  performance  of  the  tracking 
instrument,  or  the  reliability  of  the  recorded 
field  information  or  the  propagation  of  error  in 
the  statistical  techniques  used  in  analysis  and 
development  of  estimate  measures. 

In  this  paper,  two  particular  factors  are 
investigated  in  order  to  explain  the  cause  of 
some  aberrations.  Moreover,  a theoretical  basis 
is  established  which  1)  identifies  and  explains 
such  causes,  2)  can  be  used  to  search  for  pos- 
sible means  to  eliminate  such  aberrations,  and 
3)  can  be  exploited  to  determine  rationally  the 
minimum  achievable  GDOP  to  be  used  as  a standard 
basis  for  comparative  purposes. 

The  first  set  of  factors  relates  to  the 
relative  geometric  siting  of  the  instruments  as 
well  as  the  type  of  measurements  recorded.  It 
is  shown  for  example,  that  the  utilization  of  a 
radar  system  in  target  tracking  can  yield  a data 
basis  which  cannot  be  meaningfully  aggregated  to 
form  a statistical  input  basis  for  target  esti- 
mation. In  other  words  large  GDOP  values  can 
result  from  the  possible  minimal  achievable  value 
The  second  set  of  factors  relates  to  the  selec- 
tion of  a frame  of  reference  for  azimuthal  mea- 
surement. Although,  in  practice  an  established 
frame  of  reference  is  often  used,  it  is  shown 
that  coordinate  selection  can  play  for  certain 
types  of  measurements,  particularly  those  involv- 
ing photo-optical  devices  and  radars,  a signif- 
icant role  in  multiplying  the  GDOP  by  several- 
fold from  the  rationally  established  minimum 
basis. 


I . INTRODUCTION 


In  a previous  publication  [1],  a methodology 
was  developed  to  determine  and  measure  the  error 
contribution  on  estimated  Cartesian  coordinates 
of  a target  tracked  by  a complex  of  instrumentation 
systems.  It  is  known  that,  in  general,  such 
factors  as  target  position,  site  location  of 
instruments,  type  of  instruments,  type  of  measur- 
able observations  from  the  instruments  and  accu- 
racy of  Instruments  affect  the  error  quantity. 

If  (S,  jf,  z)  refer  to  the  position  of  the  esti- 
mated target,  the  error  measured  as  GDOP  is  given 


GDOP  = /var[*j  + var[J]  + var[z] 


THE  NORMAL  EQUATIONS  FOR  EQUI-ELEVATION 
SITING  OF  INSTRUMENTS  RELATIVE  TO  AN 
AUXILIARY  TARGET. 


This  GDOP  criterion  plays  an  extremely  important 
role  in  ascertaining  the  validity  of  direct  measur- 
ements in  data  reduction  techniques  to  derive 
statistical  estimates  of  the  target  state. 
Simultaneous  direct  measurements  taken  in  elevation 
azimuth  and  range  from  independent  instrumentation 
sources  at  a particular  Instant  of  time,  when 
aggregated  and  reduced  using  well  known  statistical 
techniques,  often  yield  sensible  estimates.  That 
is,  the  GDOP  associated  with  these  estimates  is 
significantly  consistent  with  the  empirical  GDOP 
values  obtained  from  past  computations.  However 
strong  aberrations  and  inconsistencies  in  the 


Consider  a situation  where  n instruments 
(observation  sites),  n 3,  of  the  same  type  are 
located  on  a circle  of  radius  d in  the  x-y  plane. 
The  system  of  coordinate  axis  is  so  selected  that 
the  position  of  an  auxiliary  target  Pc  (an 
approximate  target  position  whose  coordinates  are 
known)  is  located  on  the  z-axis  passing  from  the 
center  of  the  circle  and  normal  to  its  plane. 

The  altitude  of  P0  measured  from  the  x-y  plane  is 
h;  thus  the  coordinates  of  Pn  are  (xq,  yQ,  zQ)  = 
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(0,  0,  h) . All  elevation  angle  measurements 
of  the  auxiliary  target  P0  from  all  sites  are 
equal  to  E.  Let 

(x,  y,  z)  - coordinates  of  the  true  target 


(A°,  E°,  Ri) 
(A“,  e”,  r") 


(x,  y,  z)  - coordinates  of  the  true  target 
position  P, 

(A*,  E^,  R^)  - azimuth  angle,  elevation  angle  and 
range  from  site  i to  point  P 
0 (i  * 1,  2,  ....  n), 

(A^,  E^,  Rj_)  * azimuth  angle,  elevation  angle  and 
range  from  site  i to  point  P0, 

(A™,  E™,  R™)  “ actual  observed  or  measured  azimuth 
angle,  elevation  angle  and  range. 

Each  observation  site  is  assumed  to  operate 

independently  of  any  observation  sites.  For 

convenience,  we  introduce  the  following  quantities: 

Ax  = x - x„  AA,  = A,  - A° 

o ill 

Ay  - y - yQ  (l)  ae1  = ex  - e°  (2) 

Az  = z - z AR  - R.  - R° 

o i i i 

where  i ■ 1,  2,  ...»  n.  Each  of  the  above  A 
quantities  is  unknown  since  they  involve  non- 
measurable  variables;  however,  the  quantities 

A'Aa  - A™  - A° 


A'Rj^  - R“  - R° 

are  known  exactly.  We  assume  that  1)  {A' A, }, 
{A'E^}  and  {A'R^}  are  sample  observations  taken 
respectively  from  the  normal  populations 
N(AAt,  o|),  N(AE1,  o|)  and  N(ARi,  o|) , 2)  {A'A^,' 
{A'E,},  and  {A'R^}  form  a sequence  of  indepen- 
dently distributed  random  variables  and  3)  these 
random  variables  do  not  depend  on  the  coordinates 
of  any  of  the  sites  from  which  the  measurements 
are  made.  The  quantities  o^,  a R,  and  oR  are 
known  characteristics  of  the  instruments. 

As  shown  in  [1],  a likelihood  function  can 
be  formulated,  and  its  maximization  will  yield  the 
following  normal  equations  for  Ax,  Ay,  A~z,  the 
estimates  of  Ax,  Ay,  Az: 


Ax.l  r , 2.U  . 

T2(~2  lsin  Ai  + 


2.0  . sin  E cos  E 


2 _i-n 


d o,  i-1 
A 


V 2i° 

• l cos  A 

i-1  1 


. d cos  E r 2 .0, 

+ 2 cos  Ai] 

Or  Li 

Ajr,  1^  . sin2  E cos2  E d^_cos^_E 

O'  1 ^ o'f 


sin  A^  cos  A® 

. Az  / sin  E cos^  E . d2  sin  E cos  E , *”n 
d2  a2  a2  i-1 

cos  A0 
i 

1.  1 1rn  . .0  ...  sin  E cos  E 1"n 

’ 2 / 8ln  Ai  & Ai ~2 l 

a.  i-1  a_  i-1 

^ i“n  ® 

cos  A^  A'E^  + - [ cos  A*j*  A'RjJ  (4) 

Or  i-1 


Ax.  1 sin^  E cos2  E M d2  cos2  E^y" 

2'~  2 + ~5  + 5 ’ A 

A n -*■  4.1 


4 .0  .0 

sin  A^  cos  A^ 


. Ay  , 1 r 2 .0  . sin  E cosZ  E 
+ 2 [_2  A COS  Ai  + 2 A 


o i-1 
A 


, 2 t0  ^ d cos  E r 2 40, 

sin  + 2 L s*n 


**  3 2 i*n 

^ A z,  sin  E cos  E . d sin  E cos  Ev  ^ 

+ 2^~  2 * 2 } L 

d oE  0R  i-1 


) I Sin  At 


1,  1 r A0  A 1 A 

jt—  A cos  Ai  6 Al  ' 


\ 1’1 


sin  E cos  E. 
2 ’ 


sin  A^  A'E,  + **  \ sin  A*?  A'R  ] 

i i 2 i i 


Ax.  sin  E cos 
2i  ” 2 

d °E 


E . d2  sin  E cos  E, 1r"  .0 

- + j >/,C08  Ai 

aR  i-1 


. Ay,  sin  E cos^  E . d2  sin  E cos  E,  *>n  . .0 

+ -|( 2 + 2 > A sin  Ai 


" 4 .2,2,,  i-n 

Az.cos  E d sin  E . r . 

* 2 ' 2 + 2 ’ 4 1 


, 2 „i-n  . , „ i-n 

+ - j.  A'V 

°E  1-1  °R  11 


Note  that 

d - h cot  E (7) 

The  simultaneous  solution  of  equations  (4), 

(5)  and_  (6)_  will  yield  the  values  of  the  esti- 
mates Ax,  Ay,  and  Az.  These  estimators  are 
unbiased,  that  is,  their  corresponding  mean  is 
respectively  Ax,  Ay,  and  Az.  Further,  it  is  well 
known  that  their  joint  distribution  is  normal. 
Formally,  the  variance-covariance  matrix  can  be 
computed  once  expressions  for  Ax,  Ay  and  Az  are 
obtained.  In  particular,  explicit  formulas  for 
the  variance  can  be  derived.  For  notational 
purpose  we_  let  o2  = var[Ax],  o2  - var[Ay]  and 
a ^ - var[A~z].  Then 

2 2 2 2 

GDOP  - a + oz  + (8) 

x y z ' ’ 

In  what  follows,  we  shall  analyze  and  investigate 
the  characteristics  of  GDOP  for  two  special 
classes  of  problems: 

1.  In  the  first  case,  we  assume  that  n, 
n >_  3,  identical  instruments  having  angle  and 
range  measuring  capabilities  (e.g.  radar),  are 
sited  equidistantly  from  each  other  and  are 
recording  simultaneously  the  position  of  the  tar- 
get. 


. J. 
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2.  In  Che  second  case,  we  assume  that  n - 4 
identical  instruments  having  only  angle  measuring 
capabilities  (e.g.  cinetheodolltes,  telescopes) 
are  located  at  the  vertices  of  a rectangle  and 
are  again  recording  simultaneously  the  azimuth 
and  elevation  of  the  target. 

For  each  of  these  problems  we  shall  study 
in  particular  the  impact  of  the  elevation  angle  E 
upon  GDOP.  In  addition,  for  the  second  problem 
the  relation  between  the  selection  of  the  frame 
of  reference  upon  GDOP  will  be  investigated. 

HI.  EFFECT  ON  GDOP  WHEN  TRACKING  IN  AZIMUTH, 
ELEVATION  AND  RANGE  WITH  n >.  3 IDENTICAL 
INSTRUMENTS 

The  configuration  of  the  Instruments  siting 
is  as  shown  in  Figure  1 


i-n 


l sin(^  - a)  - 0 - l cos(^i  - a) 
i-1  " i-1  n 

1rn  , 2,2iri  . n 4?*  2,2iri  . 
I sin  (— «)  ■ - * } cos  ( - a) 

i-1  n 1 i-1  " 


i-n 


,2ni 


(10) 


i-n 
I sin( 
i-1 


2*1  • . ,2*1  . 

a)  cos( a)  » 0 

n n 


for  all  n >_  3 and  all  a 0 

2 
a 

y 


2 2 

Due  to  symmetry,  a - a . The  expression  for 
(GDOP)2  is  [1]:  x y 


(GDOP) 


A cos  E 


-|sin2E+|sin^Ecos2E  +-.cos^E 
°A  °E  aR 


FIGURE  1: r SITE  LOCATIONS  ON  THE  x-y  PLANE,  AND 
AZIMUTH  MEASUREMENT  UNDER  SYMMETRICAL 
SITING  OF  INSTRUMENTS 


We  assume  here  that  all  sites  are  located 
equldistantly  from  each  other;  then,  in  general 

A®  - 2iri/n,  1 - 1,  2 n (9) 

The  azimuth  angles  are  thus  selected  as  integral 
multiples  of  a least  value  generating  a symmetrical 
configuration  of  the  site  geometry.  As  pointed 
out  in  [1],  for  such  a geometry,  the  various  inter- 
playing terms  in  equation  (A),  (5)  and  (6)  affected 
by  the  azimuth  angles,  remain  invariant  under  a 
rotation  of  the  x-y  coordinates.  More  explicitly, 
if  a is  the  amount  of  counterclockwise  rotation 
Imparted  to  the  x-y  coordinates,  then 


1 4 2 

2 sin 


E cos  E + 


“h2" 


in2  E 


(ID 


where  E,  0 < E < -r,  is  the  common  elevation  angle 
of  the  Instruments.  It  is  evident  that  E is  a 
decision  variable;  the  Instruments  could  be  sited 
in  such  a way  that  a preassigned  value  of  E can 
be  selected,  subject  to  any  constraints  dictated 
by  physical  or  other  considerations. 

It  was  shown  in  [1]  that  when  the  instru- 
ments have  1)  angle  measuring  capability  only 
(oR  * “)  or  2)  range  measuring  capability  only 


(°a 


“) , optimum  elevation  angles  can 


be  obtained  for  each  case  which  minimizes  GDOP. 

In  the  first  case,  E(opt)  - 5A°51’,  while  for 
the  second  case  E(opt)  - 35°16'.  Neither  of 
these  values  depend  on  the  target  altitude  h 
or  the  instrument  accuracies  oA,  or,  oR. 

On  the  other  hand,  it  is  easily  seen  from 
(11)  that  when  the  Instruments  have  simultaneous 
angle  and  range  measuring  capabilities  (o a,  eg, 
o_  < ”)  the  optimum  value  of  E,  if  it  exists, 
will  depend  on  the  target  altitude  and  the  instru- 
ments accuracies. 

We  consider  now  a common  situation  when 
°A  “ °g  an<*  define  the  quantity  a - o^h/oR. 
Expression  (11)  can  be  written  as: 


(GDOP) 


°R 

Tf 


A cos  E 


2 

sin  E 


A 2 
sin  E cos  E 


+ cos  E 


sin  E cos  E 


(12) 


sln2E 


For  a - AS,  Figure  2 shows  the  plot  of  (GDOP)  as 
a function  of  E.  The  behavior  of  o2  is  regular: 
as  E varies  from  0 to  i,  decreases  from 
infinity  to  1.  The  behavior  of  cj£  and  a?  are 
quite  unexpected:  both  reach  a maximum  value  for 
an  elevation  angle  of  about  55°. 

The  plot  of  the  (GDOP)2  function  for  various 
values  of  a2  is  exhibited  in  Figure  3.  Patcerns 
of  instability  as  reflected  by  high  GDOP  values 
can  be  distinctly  noticed  on  the  graph.  Also, 
practical  considerations  often  limit  the  instru- 
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FIGURE  2:  (GDOP)2  AS  A FUNCTION  OF  THE  ELEVATION 
ANGLE  FOR  a =>  aAh/oR  - «'To7  MEASURE- 
MENTS TAKEN  IN  AZIMUTH,  ELEVATION 
AND  RANGE. 


FIGURE  3:  (GDOP)2  AS  A FUNCTION  OF  THE  ELEVATION 
ANGLE  E FOR  DIFFERENT  VALUES  OF  a. 
MEASUREMENTS  TAKEN  IN  AZIMUTH, 
ELEVATION  AND  RANGE 


ment  measuring  capability  to  elevation  angles  of 
less  than  80°.  Thus,  as  a general  rule  of  thumb, 
elevation  angles  E should  be  selected  from  the 
range  30°-50°  In  order  that  the  estimates  of 
the  Cartesian  coordinates  of  the  target  yield 
GDOP  values  which  are  close  to  the  minimum 
achievable  level. 


IV.  EFFECT  ON  GDOP  WHEN  TRACKING  IN  AZIMUTH 
AND  ELEVATION  WITH  n - 4 IDENTICAL 
INSTRUMENTS . 

The  configuration  of  the  Instruments  siting 
is  as  shown  in  Figure  A 


FIGURE  4:  SITE  LOCATION  ON  THE  x-y  PLANE  FOR 

THE  4-STATION  PROBLEM  (6  - angle 
at  the  center) . 


In  general  the  orientation  of  the  x-y  coordi- 
nates will  affect  the  estimates  of  Cartesian 
coordinates  of  the  target.  Thus,  In  (4),  (5)  and 

i-n  n i-n  i-n  . n 

(6)  the  terms  £ sin  A1:,  £ cos  aV,  £ sin  a”, 
i-1  i-1  1 i-1 

i-n  2 Q i-n  Q Q 

£ cos  A.  and  £ sin  A,  cos  A.  will  depend  not 

i-1  1 i-1  1 X 

only  on  the  relative  positions  of  each  site  1 

with  respect  to  each  other,  but  also,  on  their 
relative  position  with  respect  to  the  x-y 
coordinates.  In  turn,  the  selection  of  the  x-y 
coordinate  system  affects  GDOP.  It  has  been 
verified  [2],  that  if 
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1)  (£  cos  2 A ) - ( £ sin  2 A ) - constant 

i-1  1 i-1  1 

i“n  i-n  - 

and  £ sin  A.  - 0 - £ cos  A^  (13) 

i-1  i-1  1 

2)  or  if  readings  are  taken  in  elevation  and  range 
only,  i.e.  azimuth  measurements  are  absent  (ojr, 

Op  < “>,  o^  - -)  then  the  GDOP  will  be  Independent 
of  the  positioning  of  the  x-y  coordinates.  Neither 
of  these  conditions  are  satisfied  in  our  present 
example.  We  are  Interested  in  determining  the 
amount  of  rotation  to  be  imparted  to  the  x-y  axes 
for  minimum  and  maximum  GDOP  values.  If  6 is  the 
angle  at  the  center,  then,  the  azimuth  angles  A~, 
Aj,  Aj,  and  A^  are  given  by 


2 2 2 i,  2 

_ 1 . sin  Ecos  E , d cos  E 
M ,+  j + 2 

A QE  CR 

2 2 2 2 
_ 1 . sin  Ecos  E . d cos  E 

N + 2 + 2 

°A  °E  °R 

4 2 2 

„ _ cos  E , d sin  E 

n2  7~ 

°E  °R 

Then,  it  can  be  shown  that  [2]: 


'M(M2  - 2MNcqs29  +N2cos28) 

(M2  - N2cos26)2  4K 


M2  - N2cos26  4K 


The  plot  of  the  above  function  as  the  elevation 
angle  varies  from  0 to  ir/2,  for  various  values 
of  0 is  shown  in  Figure  5 


A.  * A + 0 + it 
A 


It  can  be  shown  that  the  amount  of  rotation  a to 
be  imparted  to  the  x-y  axes  for  least  GDOP  is 
given  by  (0  i ir/2) 

0 

l sin  2 A 
i-1 

tan  2 a - -j — ; 

£ cos  2 A.T 
i-1  1 

_ 4 sin(2A  +0)  cos  0 . 


A cos(2A  +0)  cos  0 


tan(2A  + 0) 


The  complete  details  of  the  analytic  arguments  may 
be  found  in  [2]  and  will  not  be  repeated  here. 

The  solution  of  (1A)  yields 

O If  TT 

“ " A + 2 + 2 ’ Wher"  k “ °>  41*  ±2* 

It  is  easy  to  verify  that  this  corresponds  to 
selecting  the  coordinate  axes  parallel  to  the 
sides  of  the  rectangle.  This  rotation  is  valid 
so  long  as  0 <_  0 < ir/2  and  ir/2  < 0 <_  x . When 
6 - ir/2,  the  expression  for  solving  a is  indeter- 
minate. Actually,  this  corresponds  to  the  special 
case  when  the  rectangle  shapes  into  a square,  and 
as  previously  noted  from  expression  (13),  the 
GDOP  for  that  case  is  independent  of  the  orien- 
tation of  the  x-y  coordinate  system.  If  such  an 
optimum  orientation  is  selected,  then 

A - - f + Y , k - 0,  ±1,  ±2,  ... 

In  a simlliar  fashion,  one  may  show  that  the  worst 
instrument  siting  which  maximizes  GDOP  occurs  when 

A - f - | + -y,  k - 0,  ±1,  ±2,  ... 


0 10*  20*  10*  *0*  $0*  *0*  70*  80*  90* 


FIGURE  5:  (GDOP)  /(GDOP)  AS  A FUNCTION 
max  min 

OF  THE  ELEVATION  ANGLE  E FOR  THE 
4-STATION  PROBLEM.  MEASUREMENTS 
TAKEN  IN  AZIMUTH  AND  ELEVATION 
(0  - angle  at  the  center). 

It  is  Interesting  to  note  that  for  low  values  of 
0 and  low  elevation  angles  E,  the  GDOP  value  is 
quite  sensitive  to  the  selection  of  the  x-y 
coordinate  axes.  This  sensitivity  disappears 
for  all  practical  purposes  when  the  value  of  the 
angle  at  the  center  0 exceeds  40°. 
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V.  CONCLUSIONS 

The  effect  on  estimation  of  the  Cartesian 
coordinate  position  of  a target  by  a system  where 
data  measurements  are  gathered  simultaneously  in 
azimuth,  elevation  and  range  (e.g.  radars)  is 
studied.  It  is  assumed  that  the  number  of 
measurements  exceeds  the  minimum  number  required 
thus  leading  to  an  overdetermined  system  of 
equations.  The  method  of  maximum  likelihood 
estimation  is  used  to  determine  analytically  the 
net  contribution  to  GDOP  of  the  various  inter- 
playing variables.  It  is  shown  that  low  GDOP 
values  result  when  the  measuring  instruments  are 
positioned  so  that  their  approximate  elevation 
angle  to  the  target  is  in  the  range  30° -50°. 

In  another  example,  it  is  shown  that  for 
measurements  involving  azimuthal  data  and  for 
certain  types  of  instrument  siting,  coordinate 
selection  affect  the  GDOP  values  and  that  a 
minimum  GDOP  can  be  achieved  by  properly  select- 
ing the  orientation  of  the  coordinate  frame  of 
references . 
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